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ABSTRACT 

This paper presents a completely analytic treatment of cosmological fluctua- 
tions whose wavelength is small enough to come within the horizon well before 
the energy densities of matter and radiation become equal. This analysis yields a 
simple formula for the conventional transfer function T{k) at large wave number 
fc, which agrees very well with computer calculations of T{k). It also yields an 
explicit formula for the microwave background multipole coefficient Ci at very 
large I. 

Subject headings: cosmic microwave background — dark matter — early universe 

1. Introduction 

The transfer function gives the wave length dependence of the growth of perturbations 
in the cold dark matter density from early times to near the present. As such, it plays a 
central role in theoretical studies of cosmological structure formation, and it also enters in 
the calculation of the microwave background anisotropics of large multipole number. For 
general wave length the transfer function can only be calculated numerically. This paper will 
present a purely analytic solution of the equations governing the evolution of perturbations 
in the early universe in the case of small wave length,^ which yields a simple formula for the 
transfer function in this case, including the numerical parameters appearing in this formula. 



^In speaking of small wavelengths, it is nevertheless assumed that the wavelength is large enough so that 
the fluctuations are far outside the horizon during the era of electron-positron annihilation, and large enough 
so that viscosity and heat conduction are negligible until close to the time of recombination, as is the case 
for all fluctuations of physical interest. 
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The most closely related previous work seems to be that of Hu & Sugiyama (1996). 
In contrast with their work, the present paper provides the justification for a crucial step 
in calculating fiuctuations in the dark matter density (sec footnote 4 below); it is entirely 
analytic, even in following perturbations through the era of horizon crossing and in analyzing 
the case of infinite wavelength (which is needed to normalize the transfer function); and 
explicit formulas are given for the numerical parameters in the transfer function at small 
wavelength and for the cosmic microwave background multipole coefficient Q at large I. 



Wc consider the contents of the universe to consist of radiation plus cold dark mat- 
ter plus baryons (electrons and nuclei). Wc include neutrinos in the radiation, neglect- 
ing the anisotropic part of their energy-momentum tensor, which makes possible a purely 
analytic treatment. As usual, the cold dark matter is taken to have zero pressure and 
only gravitational interactions. For simplicity at first we will assume local thermal equi- 
librium, so that the fractional changes in the baryon and radiation densities are related by 
^Pb/ Pb — ^Spr/Apr = Sr, which is a good approximation until late in the matter-dominated 
era, and we will ignore the effects of curvature and a cosmological constant, which are neg- 
ligible until near the present. Later these effects and departures from equilibrium will be 
taken into account where they are relevant. 

The evolution of compressional cosmological perturbations under these assumptions are 
governed by the equations:^ 



Here q is the co-moving wave number; a{t) is the Robertson-Walker scale factor; UR{t) 
is the radiation velocity potential; 60 is the fractional change Spu/po in the dark matter 
density pu; and dots indicate ordinary time derivatives. We are using a synchronous gauge, 
with vanishing time-time and time-space components of the metric perturbation Sgf^i,, and 
with the remaining gauge freedom removed by requiring that the cold dark matter velocity 
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Generalities 




(2) 




(3) 



^These equations are a simple extension of Eqs. (15.10.50), (15.10.51), and (15.10.53) of Weinberg (1972) 
to the multi-fluid system considered here. 
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vanishes. In this gauge, all effects of gravitational perturbations for compressional modes 
are contained in the field tlj{t) = d{5gkk{t) /2a?{t)) /dt. 

For general wave numbers these equations are too complicated to be solved analytically. 
However, for large q we can divide the evolution of the fiuctuations into two overlapping 
eras, in each of which there are approximations available that allow an analytic solution. 



3. Radiation Dominated Era 

First, consider an era so early that po and pb are much less than p^, though the 
wavelength may be inside or outside the horizon. Here a oc y/i and t'^ — 3/327rG'pi?, and by 
eliminating Ur we obtain a pair of coupled equations for Sr and ip: 

The linear combination of the three independent solutions that grows most rapidly for small 
time is^ 

2N [2 f 2\ 2\ 

^R^-Qi \e^^''^- "^""^ ~¥) ' 

^ = Ar(^^sin^+lcos^-A_l^ (6) 

where N is an unknown function of q that is presumably fixed during the era of infiation; 
9 = C\/t\ and C is the constant C = [2q^/i/ V3a]t^o. Also, Eq.(2) gives 

2N f'^^ / 4 4 4 2 \ 

Note that the fractional perturbations 6d and 6r are both of order ip/q^, justifying the 
neglect of the matter term in Eq. (1) when pR 3> Pd- 

For convenience later, it is useful to normalize the Robertson-Walker scale factor so 
that a = 1 at the time tsQ when the matter density pm = Pd + Pb and the radiation density 
Pr have a common value peq- Then at early times we have a — > (327rGpEQ/3)^/*^v^, and 
so C = (g/\/3)(27rGpEQ/3)~^^^. Also, q is now defined as the physical wave number q/a at 
t — iEQ- 



^Aside from normalization, this solution is equivalent to that given for the Newtonian potential in a 
different gauge in Eq. (48) of Bashinsky & Bertschinger (2002). 
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4. Deep Inside the Horizon 

Following this is an era in which the dark matter density may not be negligible, but 
the wavelength is well within the horizon. With the physical wave number q/a much greater 
than the expansion rate, there are two kinds of normal mode, that can be calculated using 
two different methods of approximation. 

The first are the "fast" modes, for which d/dt acting on perturbations gives factors 
of order q/a. Inspection of Eqs. (l)-(3) shows that there is a solution with 5r = 0{qil!), 
Ur = 0{ijj), and 5d = 0{ip/q), so that we can neglect the term on the right-hand side of 
Eq. (2), and even for po > Pr we can neglect the dark matter term on the right-hand side 
of Eq. (1). Ehminating Ur then gives an equation for 6r alone: 

where R = 3pb/4:P^. This has the well-known WKB solutions (Peebles & Yu, 1970) 



ij^ = (1 + fl)-"Vxp I ±, I -7=^^ I . (9) 



which would be exact for vanishing R. 

Then there are "slow" modes, for which d/dt acting on perturbations gives factors 
of order 1/t. Inspection of Eqs. (2)-(3) shows that in this case there is a solution with 
5d = 0{il}), i/j ~ q^UR. and 6r = 0{Ur) = 0{ijj/q'^). It follows that even for p^ < Pr we 
can neglect the radiation term on the right-hand side of Eq. (1), so that after eliminating 
the field ip we have 



It is convenient to convert the independent variable from t to a, using the Friedmann equation 

-^ = ^-{pM+ Pr) = ^-^(« ) (11) 

so that Eq. (10) reads^ 

^d'^So 3a\ ddo 3,^ ^, ^ „ 



"^Hu and Sugiyama (1996) pointed out that this equation leads to a transfer function with the asymptotic 
form Ink/k^, but it has not been clear why it is legitimate in deriving Eq. (12) to neglect fluctuations in 
the radiation energy density as a contribution to the source of the gravitational field during the radiation- 
dominated era. Eq. (12) was first derived by Meszaros (1974), who simply ignored fluctuations in the 



- 5 - 



where P = Pb/pm — ^b/^m- The independent solutions of Eq. (12) for /3 = were given 
by Meszaros (1974) and Groth & Peebles (1975): 

Hu and Sugiyama (1996) have given the solutions for general f3 in terms of hypergeometric 
functions, but the necessity of matching these solutions to those that apply after recombina- 
tion leads to an extremely complicated formula for the transfer function, which obscures the 
dependence of the result on the baryon density. Here we will assume that P is small, though 
not entirely negligible, and work with solutions valid only to first order in (3. The first-order 
solutions of Eq. (12) with the same behavior for a <^ 1 as the zeroth order solutions /i and 
/2 are: 

<5g'^)(a) = - f £ [hia)Mb) - h{a)Mb)] ^-^^ . (14) 

By applying Eqs. (2) and (3), we can find the fast mode solutions for and (5^ from 
Eq. (9) and the slow mode solutions for U^'"^^ and (5^'^^ from Eq. (14). These four modes 
a complete set of solutions of the fourth-order system of equations (l)-(3) up to the time 
of recombination for g/a ^ a/a. The physical solution is a linear combination of these four 
modes, to be found by matching their behavior for a 1 to that found in Section 3. 



radiation density. Groth and Peebles (1975) neglected fluctuations in the radiation density on the grounds 
that the wavelength is much less than the Jeans length, which for radiation is the horizon, but the relevance 
of the Jeans length in an expanding universe containing both radiation and dark matter is not clear. In 
their Appendix B, Hu and Sugiyama (1996) neglected the contributions of perturbations in the dark matter 
density to the gravitational field at early times, and showed that then the contribution of perturbations in the 
radiation density to the gravitational field are also negligible. But this does not justify the use of equation 
(12). For this, it is necessary to show that perturbations in the radiation density make negligible contributions 
to the gravitational field when the contributions of the dark matter perturbation are not negligible, as is the 
case late in the radiation dominated era and during the cross-over from radiation to matter dominance. Liddle 
and Lyth (2000) on p. 107 attempted to explain the neglect of perturbations in the radiation energy density 
in Eq. (12) by claiming that Silk damping makes these perturbations decay away. This is incorrect. For 
wavelengths of physical interest Silk damping is negligible during the radiation-dominated era and through 
the time of radiation-matter equality. (This has been acknowledged by Liddle and Lyth in an erratum: 
star-www.cpes.sussex.ac.uk/ andrewl/infbook/errata.html.) The neglect of perturbations in the radiation 
density in Eq. (12) is explained by counting powers of 1/q as done here, and it applies only to the slow 
mode part of the solution; in the fast mode it is the perturbations in the dark matter density that become 
negligible for small wavelength. 
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5. Matching 



Fortunately, for small wavelength there is an overlap in the two eras in which we have 
found solutions for 5di etc., satisfying both conditions q/a ^ a/ a and pM <^ Pr- In this 
period Cy^ ^ 1, and Eq. (5) gives the oscillating part of the fractional perturbation in the 
radiation density as 5r = — {2N / C'^) cos Cy/t, which for pM <^ Pr fits smoothly with the 
linear combination of the fast solutions (9) for q/a ^ a /a: 



rfast 

Or — 



2N 



cos 



qdt 



2N^6tiGpeq 

:i + i?)V4C2-"\^yo y/S{l + R)aj ~ (l + i?)VV 
from which we also find, to leading order in 1/q, 



cos 



qdt 



y/3{l + R)c 



fast 



3iV»(l + /?)-'-/H2 + /?)v/2^G'pEQ 
qH^ 



sm 



qdt 



y/3{l + R)a 



.fast _ ^Na\l + Rf/\2 + R)./6^^G^ [ qdt 



+ R)a 



and 



^fast _ 2iVyM^ 



sm 



qd.t 



(15) 



(16) 



(17) 



(18) 



To find the coefficients in the slow modes, we note that the limit of Eq. (7) for Cy^ 3> 1 



IS 



AN 

C2" 



aq 



+ 7 + ln' 

2 s/SnGpEQ 



(19) 



where 7 = 0.5772 ... is the Euler constant. For a <S 1, the solutions (14) become 



^ -ln(a/4) - 3 . 



(20) 



The linear combination of these solutions that fits smoothly with Eq. (19) is then 



rslow ■ 

On = 



'D 



AN 



7 , /2V3C2 

-- + 7 + In 

2 y q 



Od 



ANy/QnGpEQ f 



-1 + 7 + lnf ] 
2 \,/h^G^) 



S^^^ - S 



(2) 



(21) 
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The slow part of the velocity potential and radiation density are given by Eqs. (2), (3), and 
(21) as 

fj^low ^ ^slow^^2 ^ _(jgow^^2 ^22) 

Sf'- = -3a'±[a{l + R)Ut-] . (23) 

The full solution up to the time of recombination is given hy So = + and likewise 
for Ur and Sr. 

Eq. (17) shows that the fast part of is smaller than the slow part (21) by a factor of 
order l/qH"^ , so that for small wavelengths the full perturbation to the dark matter density is 
given by Eq. (21) from the time that q/a becomes much greater than a/ a, and even after the 
energy densities of matter and radiation become comparable, up to the time of recombination. 
But this is not true of the radiation perturbations. Comparison of Eqs. (22) and (23) with 
(18) and (15) shows that for large q the perturbations to the radiation velocity potential and 
density are dominated by the fast mode, by one and two factors of q, respectively. 



6. The Transfer Function 

The transfer function T{k) is properly defined as the growth of the total matter density 
perturbation for a given present physical wave number k = q/aito) = q{l + 2;eq), from early 
in the radiation-dominated era to late in the matter dominated era, relative to the growth 
that occurs in the same time interval for zero wave number. We must therefore now project 
the solution we have found for the density perturbations forward into the era following the 
time of recombination. In this era the baryonic perturbation is no longer suppressed by 
radiation pressure, and so it follows the same equation as the dark matter perturbation: 

where bu = ^Puj Pm = (1 ~ iS)bD + i^^B- This does not mean that bs and bo are equal, for 
they satisfy different initial conditions at recombination. But from a linear combination of 
these equations for bn and 5^ we find that 

, .d'^SM . A , 3a\ dbM 3 

This has the solutions (13). To find the correct linear combination of these solutions, we 
note that Sb and Sb vanish to leading order in 1/q at recombination, so Sm and and its 
first derivative at recombination must respectively equal {1 — (3)Sd and its first derivative. 
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The total matter density perturbation after recombination is the linear combination of the 
solutions (13) that matches in this way with the solution (21): 



i„(„) = '^^^y^ (^./.(a) + A,Ma)) , (26) 



where 



7 

+ 7 + ln 



2 y^/27rGpEQ J 

A2(q) = -l + (3-(3Iu-PIn 



(l-/3-/3Xi2)-/3X22, (27) 



7 

h 7 + In , , 

2 \^J2^:GpEq) 



(28) 



Ma)fjia) da 
VI + a 

Near the present, where a ^ 1, the matter density fluctuation goes to 



X,- ^ I f . (29) 

^ Jo 



l%A^Na 27rGpEQ 

^ g2 V § ■ (30) 

We can flnd the behavior of Sm early in the matter-dominated era by taking C\/t <S 1 in 
Eqs. (5) and (7): 

Sm Nt/2 . (31) 

Eqs. (30) and (31) must be compared with the growth of Sm for q — 0. In this case 
Eq. (2) gives Sr = Sd = —ip, so Eq. (1) becomes 

d I 2 dSM \ , ^ o /I 



The solution of Eqs. (32) and (11) that has the same behavior for a — > as Eq. (31) is 

N 



JM 



5a2 Y 27rGpEQ 
This has an asymptotic behavior for a ^ 1: 

Na 



(iQ + 8a - 2a^ + - leVTTa) (33) 



JM 



5 y 27rGpEQ 



(34) 



The transfer function T then has an asymptotic behavior for large wave number given by 
the ratio of Eqs. (30) and (34): 

T^'^AM), (35) 



-9- 



with Ai{q) given by Eq. (27). At late times the growth of Sm may be affected by a cosmo- 
logical constant or spatial curvature, but these effects are independent of wave number, and 
therefore cancel in the transfer function. 



For <^ pd we can neglect /3, so Eqs. (35) and (27) give a transfer function 



T 



q2 



7 , / 2q 

h 7 + In , 

2 \,/2^^G^ 



(36) 



This can be put in more famihar terms by using the relations peq = {^H^VIm /9>t^G) (1+^eq)^, 
q — A;(l + -Zeq), and 1 + 2;eq = ^m/^r, which give the transfer function in terms of the 
present wave number k: 



-- + 7 + ln 



ln(2.40Q) 
(4.07Q)2 



(37) 



where Q = A:(Mpc^"'^)/r2jv//i^, and in the final expression wc use ClRh"^ = 4. 15 x 10^^. This may 
be compared with the BBKS numerical fit (Bardeen et al. 1986) to computer calculations 
of the transfer function: 



T{k) ~ 



ln(l + 2.34g) 
2.34g 



[1 + 3.89g + (16.ig)' + (5.46g)3 + (6.71Q)^] 



-1/4 



(38) 



This goes to ln(2.34(5)/(3.96(5)^ for large Q, in very good agreement with Eq. (37). Our 
simple calculation thus accounts not only for the form of the transfer function for large wave 
numbers, but also for its numerical parameters. 

(Though it is not relevant to the present work, it may be noted that the BBKS formula 
cannot be taken seriously for small values of Q, for it has unphysical terms that are linear 
in Q at Q — > 0. Analyticity in the three- vector k requires that in this limit T{k) should be 
a power series in /c^, or equivalently in Q^.) 

To assess the effect of a non-zero baryon number, we note that, to first order in (3 — 
^b/^m, the general formula (35) may be put in the form 



T 



In (^2.40g(l + /3X22; 
4.07q(i + /3(1+2:i2)/2) 



(39) 



We need values for the integrals Xij defined by Eq. (29). The upper limit on the integrals 
(29) is = {1 + 2;eq)/(1 + -Zr). The redshift Zji at recombination has only a very weak 
dependence on cosmological parameters, and will be taken here to have the fixed value 



zr — 1100. The redshift z-^ci at matter-radiation equality is given by 1 + zeq = ^.m/^r, so 
taking QrH'^ — 4.15 x 10~^, we have qr — 21.9 Qm^^- The integral lu is given by 



Ii2{aR) = — -22aii - 18a^ + 4 In 




( 



vi + OR-iy 



The integral I22 is given by a lengthy expression involving Spcncc functions, but it converges 
so rapidly for likely values of gr that for practical purposes we can use the value for gr 
infinite: 



For instance, for gr = 4.38 (corresponding to Vtuh^ = 0.2), we have T22 = 0.9470. 

Eq. (39) agrees very well for large k with the numerical results of Holtzman (1989). For 
flMh^ = 0.2 (the most plausible of the values considered by Holtzman) Eq. (39) becomes 



with k in Mpc~^. Table 1 compares the results given by this formula with the numerical 
results given by Holtzman (1989) for {3 = Ub/^m equal to 0.01 and 0.1, and for various values 
of k. As can be seen, for these parameters the asymptotic formula (39) gives a pretty good 
approximation to the numerically calculated results for k > 0.5Mpc~^, and the numerically 
calculated results converge rapidly to Eq. (39) for larger values of k. But Holtzman warns 
that his result should not be used for k > 3.09 Mpc~^, while the results obtained from 
Eq. (39) presumably become increasingly more accurate for larger values of k 

We see that for any plausible value of O^uh^, the effect of a small baryon density on 
the argument of the logarithm is to replace the parameter Qm^^ in the definition of Q with 

^^^'^(1 — .95/9), while for VtMh"^ in the range of 0.12 to 0.2, the effect of a small baryon 
density on the denominator of T{k) is to replace the parameter VtMh"^ in the definition of 
Q with nufi^il - CP), with C = (1 + ^22) in the range of 1.24 to 1.38. These results 
throw some light on a series of attempts to correct the transfer function for the effects of 
baryon density by re-scaling the definition of Q (usually called q) in the BBKS formula 
(38). Various authors attempted to correct for the baryon density by replacing flMh^ in 
the denominator of Q with a factor Qm^^ exp(— 2Qb) (Peacock & Dodds 1994), or with 
Qm^^ exp(— — V^KQe/^m) (Sugiyama 1995) or with Qm^^ exp(— — Qb/^m) (Liddle 
et al. 1996). (For a more detailed study of the effects of a finite baryon-to-dark matter ratio 
on the transfer function, see Eisenstein & Hu (1998).) Of course there is no reason why 
the baryon density should enter only in the definition of Q, and Eq. (39) shows that it does 



^22(00) = 27rV5 - 3 = 0.947842 . 




- 11 - 



not. But even without any detailed calculations, it is evident that these correction factors 
are physically impossible. The transfer function is defined with no reference to the present 
moment, except that it is conventionally written as a function of the present wave number k. 
It depends on ^Mh"^ and QrH'^, which enter in the formulas for k/q and Peqj and it can (and 
does) have an additional dependence on the constant ratio of the energy densities of baryons 
and all matter, which is equal to Qb/^m, but there is no way that it can depend separately 
on or flM or h. What we have found here is that for large wave number, the effect of a 
small baryon density can be crudely taken into account by replacing the parameter Qm^^ in 
the definition of Q with Qm^^(1 — C^b/^m) — ^nh"^ exp(— CQb/Qm), with ( roughly equal 
to unity for likely values of Q.Mh'^- 



7. Microwave Background Anisotropics 



The fractional temperature fluctuation in a direction n takes the general form (apart 
from late-time effects): 



AT(n, z) 
T 



(40) 



where p{z) dz is the probability that last scattering will occur between rcdshifts z and z + dz\ 
dA^z) is the angular diameter distance to redshift z] and et is a primordial fluctuation 
amplitude, defined as proportional to -/V(k), with a coefficient to be chosen below. In the 
synchronous gauge and hydrodynamic approximation used here, and now making the further 
approximation that dark matter dominates the energy density at last scattering, the form 
factors F and G in Eq. (40) are given by (Weinberg 2001): 



eG — —aqUR + qt4>/a , 



(41) 
(42) 



where = —AirGpRSRa^/q^ is the Newtonian potential produced by dark matter density 
fiuctuations. The first and second terms in F arise from the Sachs- Wolfe effect and intrinsic 
temperature fluctuations, respectively. The form factor G arises from the Doppler effect, with 
its flrst and second terms contributed by velocities produced by pressure and gravitational 
forces, respectively. 

The conventional multipole coefficient is given in general by the familiar formula 



Ci = 16tt^ / V{k) k^ dk 
Jo 



dzp{z) [je(^kr{z)/Ho^ F{k,z) + j',(^kr{z)/Ho^ G{k,z) 



(43) 
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where V{k) is the power spectral function, defined by 

(ekek')=n^)<5'(k + k'), (44) 

and 

r{z) ^ (1 + z)dA{z)Ho = f . (45) 

To fix the normahzation of ek, we note that for small values of i (say, i < 10) the large value 
of dA makes the spherical Bessel functions in Eq. (43) oscillate rapidly except for small values 
of k, so Ci is dominated for small i by the Sachs- Wolfe term in F{k), for which Eqs. (34) 
and (41) give the z — independent small-A; behavior 

We therefore define by 

e. = - iV(q) , (47) 

5g2y67rGpEQ 

so that F(0, z) — 1. With this normahzation, a Harrison-Zel'dovich power spectral function 
V{k) = Bk-^ gives Q = S7r'^Bi{i + 1) for small £. 

For large £, the integral over A; in Eq. (43) is dominated by large wave numbers. In 
this case, the Sachs-Wolfe term in Eq. (41) receives a contribution of order 1/k^ from the 
slow mode part (21) oi Sd and of order 1/k^ from the fast mode part (17). The intrinsic 

fluctuation term in Eq. (41) receives a contribution of order l//c^ from the fast mode term 
(15) in Sr, and of order 1/A:^ from the slow mode term (23). The slow mode parts of the two 
terms in the Doppler form factor (42) cancel, leaving the contribution of the fast mode term 
(18) in Uji, which is of order We conclude from this that in the absence of dissipative 

effects, the temperature fluctuation is dominated for large k by the fast-mode part of the 
intrinsic temperature fluctuation. 

But for very large k the rapidly oscillating fast mode is killed by Silk damping (i., 
e., photon viscosity and heat conduction) and Landau damping (cancelations due to large 
changes in the phase of the fast modes over the range of redshifts at which last scattering 
may occur). As pointed out by Hu and Sugiyama(1996), for i greater than about 4,000 
the dominant contribution to Q arises from the non-oscillatory terms in the perturbations. 
These terms, which are contributed by both the Sachs-Wolfe effect and the intrinsic tem- 
perature fluctuations, can be taken from Eqs. (83) and (84) of Weinberg (2001), with the 
damped terms neglected and an extra factor T{k) supplied, because here we are dealing with 
wavelengths that come into the horizon during the radiation dominated era. This gives 

F(k, z) -3R(z) T{k) , G{k, z) ^ , (48) 
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so Eq. (43) becomes 



d = 1447r^ V{k) T^(k) k'^ dk J dz p{z) R{z) je (kr(z) /H^ 



(49) 



To do the double integral over z and A;, we use an approximation of Hu and White (1996). 
The last-scattering probability distribution p{z) is sharply peaked around a mean value zl ~ 
1. 100, while for sufficiently large I the spherical Bessel function is even more sharply peaked 
at a value £ + 1/2 of its argument. We therefore set z ai a. value where kr{z)/ Hq = i + 1/2 
everywhere but in the argument of ji, and integrate over the argument of ji with k fixed, 
after which we set k — {i + l/2)HQ/r{zL) everywhere but in the argument of p{z), and 
integrate over that argument: 



X 



dz 



r{zL) 



ji{s) ds 



(£ + 1/2)^0 
r{zL) 



RHzl) 



+ 1/2 



r'^{zL)r'{zL) 



M 



r'^{zL)(J 



V 



r{zL) 



r{zL) 



RHzl) 



where a is defined by 



P 



iz) dz 



(50) 



(51) 



20r(7 ' 

so that (7 is the standard deviation if p{z) is Gaussian. For instance, for a straight spectrum 
with V{k) oc k~^~^% Eq. (45) gives Q oc ^"^""^In^^. Unfortunately the interposition of 
foreground objects makes it unlikely that this can be measured. 

I am grateful for valuable discussions with S. Bashinsky, E. Bertschinger, R. Bond, W. 
Hu, A. R. Liddle, D. H. Lyth, H. Martel, and R Shapiro, and for help with integrals by M. 
Trott. This article is based on work supported by the National Science Foundation under 
Grant No. 0071512, and also supported by The Robert A. Welch Foundation. 
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Tabic 1. Values of the transfer function for Q^^fh"^ = 0.2 and Qb/^m = 0.01 or 0.1. 



k (Mpc-^) 




T{k)nB/nM=o.Qi' 


T{k)iig/nM=o.i'^ 


T{k)nB/nM=o.i' 


0.1 


0.161 


0.0451 


0.138 


0.0509 


0.3 


0.0398 


0.0337 


0.0328 


0.0284 


0.5 


0.0189 


0.0169 


0.0154 


0.0140 


1 


0.00640 


0.00586 


0.00517 


0.00480 


2 


0.00202 


0.00187 


0.00162 


0.00152 


3 


0.000997 


0.000938 


0.000797 


0.000762 



^Prom Holtzman (1989) 
''Prom Eq. (39) 



